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ABSTRACT 

We present an infrared study of the supernova remnant (SNR) N49 in the 
Large Magellanic Cloud with the near-infrared (NIR; 2.5 — 5 /im) spectroscopic 
observations performed by AKARL The observations were performed as a coarse 
spectral mapping to cover most of the bright region in the east, which enables us 
to compare the distribution of various line emission and to examine their corre- 
lation. We detect the 3.3 /zm aromatic feature in the remnant, which is for the 
first time to report the presence of the 3.3 fim aromatic feature related to a SNR. 
In the line maps of H2 1-0 0(3), 3.3 /im feature, and Bra, the distribution of the 
aromatic feature shows overall correlation with those of other emissions together 
with regional differences reflecting the local physical conditions. By compari- 
son with other archival imaging data at different wavelengths, the association of 
the aromatic emission to other ionic/molecular emission is clarified. We examine 
archival Spitzer IRS data of N49 and find signatures of other polycyclic aromatic 
hydrocarbon (PAH) features at 6.2, 7.7, and 11.3 /im corresponding to the 3.3 
yum aromatic feature. Based on the band ratios of PAHs, we find that PAHs in 
N49 are not only dominantly neutral but also small in size. We discuss the origin 
of the PAH emission in N49 and conclude that the emission is either from PAHs 
that have survived the shock or PAHs in the preshock gas heated by radiative 
precursor. 



Subject headings: dust, extinction — infrared: ISM: lines and bands — ISM: 
individual (N49) — ISM: supernova remnants — Magellanic Clouds 
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1. Introduction 

Over the past 20 years, various researches on polycyclic aromatic hydrocarbon (PAH) 
in interstellar medium (ISM) have been conducted both observationally and theor etically. 



PAH s are revealed as abundant, ubiquitous, and dominant in the ISM of galaxies ( iTielens 



20081 . and references therein), and one of the main resources of infrared (IR) emission from 
near- to mid- IR wavebands. There have been a number of observations to show a variety of 
PAH features at different wavelengths. Several major emission features are usually detected 
in diverse objects such as the 3.3 /xm C-H stretching band, the 6.2 and 7.7 /xm C-C stretching 
bands, and the 8.6 and 11.3 /im C-H in- and out-of-plane bending bands. 

In PAH processing, interstellar shocks are known to play a crucial role. It is largely 
accepted that PAHs could be either the remaining dust condensation nuclei that escaped 
the grain growth process in asymptotic giant branch ejecta or t he fragmentation of dust 



grain s through shattering collision in fast interstellar shocks (e.g., lTielensll2008l : I Jones et al. 



19961 ). Also, interstellar shock waves have been considered as one of the main mechanisms 
to destroy PAH molecules. In spite of plentiful observations of PAH features, however, there 
are still unanswerable questions on the role of interstellar shocks in the evolution of PAHs. 
Particularly, how supernova shocks affect PAH molecules is barely explored. Detection of 
PAH features in supernova remnants (SNRs) is unexpectedly rare, recalling that the one of 
their formation process is related to interstellar shocks. 

There have been considerable literatures dealing with the survivability of PAH molecules 
in shocked environments. In particular, since PAH emission has been categorized as one of the 



IR em ission components expected for SNRs based on the Spitzer imaging data f lReach et al. 



20061 ). substantial efforts have been made in order to search for observational evi dence of PAH 



emissi on in SNRs. By using Spitzer Infrared Spectrograph (IRS) observations, iTappe et al. 



( I2OO6I ) firstly reported the detection of 15 — 20 /im hump attributed to C-C-C bending 
modes of large PAHs (~ 4000 C-atoms) with the weakly detected 11.3 fim PAH feature in 
SNR N132D, which is located in the Large Magellanic Cloud (LMC). Based on the lack of 
PAH features at 6 — 9 fim and the large ratio of the 15 — 20 /j.m hump to the 11.3 /xm feature, 
they interpret that small PAHs are rapidly destroyed by thermal sputtering in the supernova 



blast w ave. Among the Spitzer IRS spectra of several galactic SNRs, some in iNeufeld et al. 



( 120071 ) and lHewitt et al.l (j2009[ ) also show the major PAH features with strong ionic and/or 
molecular lines, yet both authors have not mentioned the association of the features with the 
SNRs in the papers. In most cases that PAH emission is observed in SNRs, no convincing 
evidence for the PAH emission intrinsic to the SNRs has been reported. This lack of detection 
of PAH features in SNRs results from the difficulty of discrimination between PAH emission 
from SNRs and that from other back/foreground sources. 
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In this respect, the LMC staying away from the galactic disk, could be the best place 
to observe PAH emission from SNRs. Moreover, since PAH molec ules in shock regions ar e 
supposed to be able to survive only in those with dense clumps flMicelotta et all l2010al ). 
SNRs interacting with dense circumstellar or interstellar material would be good candi- 
dates for detecting PAH emission. Thus, we have targeted the SNR N49, one of the SNRs 
interac ting with ambient mol ecular clouds in the LMC, with AKARI infrared space tele- 
scope jMurakami et al.l lioOTi). N49 (SNR 0525-66.1) is a middle-aged SNR (~ 6600 yr. 



Park et al.l 



20031 ). and the surrounding dense ambient medium including complex filamen- 



tary structur es suggests that the progenitor of the SNR is a B-type star without a strong 
stellar wind flShulllll983l ). Thanks to its sufficient brightness at almost all wavel engths, this 
remnant has been studi ed by various ob servations at multi wavelengths: X-rav (IPark et al. 



20031). ultraviolet fUVl flBlair et al.ll2000r). o ptical fIVancura et al 



1992; 



Bilikova et al 



2007L 



IRjIWilliams et al.l[2006l : ISeok et al.lboosh . and radio JPickel k Milne I Il998l : IPickel et al. 



I995I ). The observational results have shown not only complex filamentary morphology in 



both optical and IR wavebands, but also a complete shell-type structure in X-ray and radio. 
By CO observations (IBanas et al.1 119971 ). N49 is suggested to be interacting with its nearby 
molecular cloud at southeast, which is thought to be responsible for mostly asymmetric 
morphology and brightness enhancement toward southeast seen at all wavelengths. In this 
paper, we present the near-IR (NIR) spectra of N49 in the unprecedented wavelength range 
(~ 2.5 — 5 yum) and report the first detection of an aromatic feature at 3.3 yum in SNRs. 



2. Observations 



We perfor med spectroscopic observations of N49 (SNR 0525-66.1) using the Infrared 
Camera (IRC; lOnaka et al.l 120071 ) on board AKARI. These observations are part of the 
several AKARI mission programs, and the details of the observations are listed in Table 
[U The data were taken with a grism in the IRC NIR channel, NG (2.5 — 5 /im), in the 
common mode of slit spectroscopy (Ns) during the post He phase (cooled by the onboard 
cryocooler). The dispersion (wavelength increment per pixel) is ^ 0.01 /im, and th e slit 
width is 5", which gives a spectral resolution of i? ~ 100 at 3 /xm (jOnaka et al.l 120071 ). As 
a coarse spectral mapping, we made fourteen pointed observations to the SNR covering a 
bright wedge-shaped feature in the east as well as relatively IR faint region at west. In order 
to obtain a background spectrum, we also carried out two independent observations toward 
neighbor regions avoiding diffuse emission from the remnant. Figure [T] shows the positions of 
the IRC NG/Ns slits superposed on the IRC N3 band (2.7 — 3.8 /xm) image that was taken 
in part of the AKARI large-scale survey of the LMC (PI: T. Onaka, llta et al.l 120081 ). The 
total integration time is 396 s for each observation. 
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The data were reduced by using th e standard IDL pack age for the IRC spectroscopy 
data reduction (toolkit version 20090211. lOhyama et aLll2007l ). General preprocesses such as 
dark subtraction and linearity correction were applied during the pipeline. The extraction 
apertures of the final spectra are determined for areas observed with more than one slit posi- 
tions, and finally five regions, named "PI" to "P5" from east to west, are chosen as marked 
in Figure [1] (white rectangles). To extract each spectrum, we integrated 21 pixels along the 
slit for source spectra (except "P5") and 17 pixels for background spectra, corresponding to 
about 31" and 25", respectively. For the "P5" spectrum, we extracted it from a shorter aper- 
ture (13 pixels corresponding to 19") in order to avoid low surface brightness regions. The 
center and the length of extraction were determined to maximize the signal-to-noize ratio 
(S/N), to avoid contamination from any point sources, and to average the resultant spectra 
in the overlapping region for the final spectra. The final spectra including the background 
spectrum are shown in Figure [21 Since the background spectrum does not show any line or 
continuum emission, we did not subtract it from the final spectra to not degrade the S/N. 



3. IR spectrum of N49 
3.1. AKARI IRC NG spectrum 

We examined the IRC NG spectra for the SNR N49, and the detected lines are given 



in Ta ble [2j The MIR emission in N49 is dominated by ionic line emission f l Williams et al. 



20061 ). However, there is no strong ionic line within this IRC NG wavelength coverage, so 
that we could not detect any ionic lines except marginal detection of [Fe ll] A 4.889 /im. The 
most prominent lines are hydrogen recombination lines such as Bra A 4.052 /xm and Br/3 A 
2.626 /im, which are seen at all spectra except spectrum PI (Fig. [2]). Since PI is extracted 
from just outside of the bright Ha region, there is no strong recombination line emission 
but only an interesting feature at 3.3 yum (see next). Although H2 molecular lines are not 
as distinct as the recombination lines, several rotational lines such as H2 1-0 0(3) and 0-0 
S(ll) lines are also clearly detected. 

For the lines detected in most spectra, we measured their intensities by fitting Gaussians 
with linear baselines. The intensities with uncertainties are given in Table [2J The uncertain- 
ties are the 1 a errors from the line fits, and the upper limits at 2 a level are given if the 
line emission is not apparent. Most lines in Table [2] are not blended with others, yet there 
are two cases (i.e. H2 1-0 0(5), Pf(5 at 3.235, 3.297 /xm and Pf/3, H2 0-0 S(9) at 4.654, 4.695 
yum) that two lines are located close to each other. For those cases, we applied two Gaussian 
fits to estimate their intensities together with either the peak wavelength or the full width 
at half maximum (FWHM) fixed if necessary. In addition, we fixed peak wavelengths or 
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FWHMs for a few relatively weak lines. 



In the case of PI to P3 spectra, we notice that the feature at 3.3 fim could largely 
originate from the well-known aromatic C-H stretching transition at 3.3 /xmrather than Pf5 
at 3.297 fim. In order to be attributable to PfS emission, its intensity is abnormally stronger 
than theoretical value (< 10% of Bra intensity in Case B) while the other Br/3 and Pf lines 
are generally consistent with the expected values (see Table [2]). Moreover, the intensities 
of the 3.3 //m feature do not vary with Bra intensities (see Fig. [2]), and the measured 
widths of the feature appear to be wider (FWHM ~ 0.04 — 0.05 /xm) than those of other 
ionic /r nolecular lines (FW HM ~ 0.03 fim on average) that corresponds to the instrumental 
width fjOnaka et al.ll2007l ). Thes e widths are consistent with t he typical width (~ 0.04 fim) 
of the 3.3 fim aromatic emission fjvan Diedenhoven et al.ll2004l ). Hence, we consider that the 
3.3 fim feature mainly arises from the aromatic hydrocarbon emission due to C-H stretching 
mode, which is the first time that we detect this aromatic feature associated with a SNR. 
Although this 3.3 fim feature does not require a polycyclic molecular structure, it is referred 
to as a 3.3 fim PAH feature en bloc hereafter. 

There are several minor PAH features near the 3 .3 fim feature such as the weak 3.4 
fim band, or the broad plateau at 3.2 — 3.6 fim (e.g., iTielend l2008l ). In the IRC spectra, 
these minor PAH features are not obvious (Fig. [2]). We could only find features at 3.4 — 3.5 
fim range from P2 and P3, but they do not look similar to each other (i.e., have different 
peak positions), and there seems no such feature at PI. We measured their intensities which 
may be regarded as an upper limit of the minor 3.4 fim PAH feature intensity consider- 
ing that other minor features and/or the underlying plateau might be included (Table [2]). 
Despite of our ambiguous detection, the presence/absence of the 3.4 fim feature could be 
an interesting issue, because this feature can give a constraint on the band carriers, which 
is still controversial. There are two interpretations for this featur e; overtones (i.e., anhar- 
monicity of hot bands; Allamandola et al. 1989 : Geballe et al.l 19891) and the C-H stretc hing 



emission from aliphatic side groups attached to PAHs (jJoblin et al.lll996l ; lTielensll2008l . and 
references therein). Excited PAH bands of C-H stretch transition can be shifted to longer 
wavelengths due to anharmonicity. In that case, the 3.4 fim intensity increases if excitation 
gets higher. Besides, smaller PAHs tend to have a higher ra tio of the 3.4 to 3.3 fim bands 
for the same internal excitation energy (IGeballe et al.lll989l ). Meanwhile, the presence of 
aliphatic bonds may delineate a balance between reactions with C, C"*^, and H producing 
aliphatic groups and photodissociation of these peripheral groups. When the UV radiation 
is strong, the dehydro genation of PAHs sev erely proceeds so that the 3.4 fim emission be- 
comes weak. However, iTielens et al.l (119941 ) suggested that ion bombardment can form the 
aliphatic bonds by hydrogenation of an amorphitized graphite surface in interstellar shocks. 
Although its detection is tentative, the possible presence of the 3.4 fim feature at P2 and 
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P3 but PI may indicate that the PAHs inside the SNR are generally smaller with a strong 
UV field or that aliphatic bonds are more abundant due to shock processing compare to 
a preshock region. Further investigation with a higher sensitivity and spectral resolution 
would be needed to understand this point more clearly. 



3.2. Brightness distribution of PAH, Ha, and H2 emission 

Although the 3.3 /im feature can be attributed to PAH emission, it is important to 
confirm whether the PAH emission is really associated with the SNR. For this purpose, 
maps of H2, PAH, and Bra emissions are constructed from all spectra covering the PAH- 
bright area in the east (Fig. [3]). The pixel size of the map is the same as the original pixel size 
of the slit (1.46" x 1.46"), and the intensity at each pixel is calculated by resampling them 
with a nearest-neighbour interpolation methoci^ and averaging the data values that fall on 
the pixel. Then, the line maps are smoothed with a three-pixel Gaussian. The resultant line 
maps reveal the distribution of 3.3 /xm PAH emission together with Bra and H2 1-0 0(3). In 
order to remove the contamination by PiS line emission, the PAH map is subtracted by the 
scaled (xO.l) Bra map. The morphology of PAH emission seemingly has association with 
those of other emission. The peak of PAH emission is found to be located at the tip of the 
bright wedge-shaped emission region, and the overall shape also shows enhanced emission 
along the wedge-shaped region. This strongly supports that the PAH emission is indeed 
related to the SNR. 

All line maps in Figure |3] show the enhanced emission at the eastern region. Previous 
observations at the other wavelengths similarly show the brightness enhancement, which is 
interpreted as the effect of interaction between the S NR and the nearby molecular cloud at 



the southeast flBanas et al.lll997l : iBilikova et al.ll2007( ) . The morphology of the PAH emission 



however, has similarities and differences in details with respect to those of H2 0(3) and Bra 
emission. The eastern boundaries of both H2 and PAH emission are indented with respect 
to the Ha boundary (e.g., a gap at the east side between the two white crosses), while the 
Bra emission shows a good agreement with it (also see Fig. [51(a) and (c)). In particular, it 
is noticeable that both H2 and PAH exceed the Ha boundary around the tip of the wedge- 
shape region, which is at the center of the interacting molecular cloud (i.e. PI and P2 in 
Fig. [T]). The two clumps of the PAH emission in the south (white crosses in Fig. [3]) could 
be seen at the coincident positions in the H2 map. Meanwhile, the peak position of the 



^The resampling uses a kernel extending a single input pixel, 
https: / /www. astromatic.net /pubsvn/software / swarp /trunk / doc / swarp.pdf 



See 
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PAH emission is not matched with that of the H2 but rather closer to that of the Bra. In 
addition, it is interesting that relatively weak PAH emission is observed where H2 emission 
is the strongest. 

Figure H] compares the average intensities of the representative emission lines in indi- 
vidual slits from PI (east) to P4 (west). Since the intensities (listed in Table [2]) are those 
averaged over the slits, they may not directly describe the spatial variation but quantita- 
tively show their trend. The peak intensity of PAH emission appears at P2 (more eastward) 
unlike others, and the emission at PI (outside the shock boundary) still has a comparable 
intensity to that of the PAH emission at P2. We could not see any tight correlations of 
the PAH emission with either H recombination or H2 molecular line emission, but it is no- 
table that the recombination line intensity rapidly increases over the boundary (P2) whereas 
the molecular line intensity is steady at the east (even inside the boundary) and becomes 
stronger at P3. 

Besides the emission intensities, we examined the variation of physical quantities such 
as the column densities, A^(H2) and the gas temperatures, Tgas with respect to the four spec- 
tra, assuming local thermodynamic equilibrium (LTE) estimated from the intensities of H2 
molecular lines. H2 line intensities can directly provide the column density, N{Ju) of each 
upper state of the transitions, Ju- An extinction correction is not applied beca use the extinc- 
tion toward N49 is negligible according to the measured E( B — V) = 0.37 fIVancura et al. 



I992I ) and the extinction curve of the "average" LMC model fjWeingartner fc Draindl200ll ). 
Performing one-temperature LTE fits to the observed H2 transitional lines with a given 
ortho-to-para ratio (OPR = 3), it is found that A^(H2) varies from ~ 1 to 5x10^^ cm"^ 
(Fig. H]), while Tgas steadily remains at ~ 2000 K. The column densities we derived are much 



lower than the typical total column densities for SNRs (~ 10^° cm~^, e.g., 
2OO9I ). This is probably because only H2 lines of high upper state transition are used for 
our estimation. In the case for SNRs, an ankle-like curve is representatively seen in the H2 



Hewitt et al. 



level populations (e.g.. IShinn et al.ll2010l ). so that H2 pure rotational lines with low upper 



states such as H2 S(0)— S(2) transitional lines mainly determine the total column density. 
Then, what we derived is likely to be the column density of a hot component (Tgas ~ 2000 
K) within the SNR. Therefore, with the current observations, we cannot explore the relation 
between total H2 and PAH. 

The morphological characteristics of PAH emission with respect to molecular and ionic 
gas would be reconfirmed and clarified by comparing to narrow filter images with a higher 
spatial resolution. As complement, we use preexist ing data, the H2 1- S(0) at 2.12 /xm taken 
with the Anglo — Austr alian Telescope {AAT) (jPickel et al.lll995l ). the S yitzer IRAC 8 



/xm (IWilliams et al.ll2006l ). and the Hubble Space Telescope {HST) Ha image (IBilikova et al. 
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20071 ). The 2.12 /im and the Ha images generally show consistent distribution with the H2 
0(3) and Bra maps in Figure [3l respectively. In Figure O we can confirm the features 
described above. The eastern boundaries of the 2.12 /xm and the PAH emission show a good 
agreement, and both emission are almost absent in the western part of the SNR where the 
optical filaments are comparably bright as seen in P5 spectrum (Fig. [2]). The peak of the 
PAH emission, however, shows better correspondence to that of Ha emission rather than 
that of 2.12 /im. It is certain that the PAH emission exists beyond the Ha boundary in east. 
Thus, the other observations of N49 similarly show the same morphological characteristics 
as our AKARI line maps. 

The characteristics of the PAH emission in N49 may be summarized as follows: (1) 
The overall morphology of the 3.3 /im PAH emission has the similarity to that of the H2 
emission in the sense that both are spatially confined. In contrast, the peak of the PAH 
emission is located rather closer to that of the Bra emission. (2) Toward the region with 
the brightest H2 emission, however, the PAH emission is relatively faint. (3) There is PAH 
emission beyond the shock front in the east i.e., beyond the east boundary in optical, and its 
peak position is near the center of the interacting molecular cloud (Fig. [T]). In addition, it is 
interesting that the distribution of PAHs shows a good agreement with the bright emission 
seen in the IRAC 8 /tm image. A noticeable feature in the 8 /im image, not covered by our 
AKARI observation, is the bright emission beyond the eastern boundary. The nature of this 
feature is described in more detail in 5 3.3. 



3.3. Spitzer IRS spectrum 

The detection of the PAH emission at 3.3 /im implies the possibility of other major 
PAH band emission such as 6.2, 7.7, and 11.3 /ir n features. We exami ned the Spitzer IRS 



archival data for PAH emission (AOR 6586112. IWilliams et al.M2006h . The IRS slits are 
centered on the brightest tip of the wedge-shape region (Fig. [5l (fe)), and partially overlap 
with the IRC slits (P2 and P3). Retrieving the IRS low-resolution (SL: 5.2 — 14.5 //m) 
Post-Basic Calibrated Data (PBCD) from the Spitzer archive, we have extracted spectra by 
using the standard SPICE package with a full aperture for an extended source (slit width: 
3.7", aperture length: 50.4" (28 pixels)). We have applied background subtraction by using 
different order spectra. One of the benefits from the background subtraction is that the 
residual of the instrumental pattern (so called "jail-bar" patterij^ seen weakly in the PBCD 
data can be efficiently removed without much loss of S/N. The final IRS spectrum is shown 



http: / /ssc. spitzer. caltech.edu/IRS/irsinstrumenthandbook/ 61 / 
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in Figure [61 It is almost the same as the one pubhshed by IWilhams et aLl ( 120061 ) , but shows 
small differences probably due to background subtraction. The strong ionic lines such as 
[Fe II] and [Ne II] at AA 5.32 /xm and 12.83 /xm, are detected together with s e veral H2 



transitional lines up to 0-0 S(2) A 12.30 fim (Fig. El (a)). As in IWilhams et all mO& \. a 
non-zero background level is seen, presumably because the background level varies. Taking 
a close look at the PAH bands in the spectra, we could identify several PAH features (Fig. |6l 
(b) — {d)). PAH emission of C-H out-of-plane bending mode at 11.3 /im is clearly detected, 
and weak C-C stretching features at 6.2 fim and 7.7 /im with ambiguous 8.6 fim emission of 
C-H in-plane bending mode can be seen. We derive the intensities of the PAH band features 
by simple summations with a linear baseline fit over the entire order. The intensities of the 
6.2, 7.7, and 11.3 fim PAH features are estimated as 1.0 ± 0.5, 2.0 ± 0.4, and 1.6 ± 0.1 in 
units of 10~^^ W m~^, respectively. The quoted errors include the calibration uncertainty 
and 1 a fluctuation in the baseline spectra. 

For the clearly detected 11.3 fim PAH feature, we compare its spatial distribution to 
those of H2 0-0 S(2) A 12.3 /im and [Ne II] A 12.8 fim along the slit (Fig. [7]). For comparison, 
the spatial profile of Ha emission extracted from the HST image is also shown. Although 
a sharp jump in the Ha intensity is not seen, a shock front might be located near the place 
where the Ha intensity rapidly increases. A dashed line in Figure [71 can be regarded as 
the SNR boundary. The emission outside the shock front could originate from ambient gas 
heated by radiative precursors, and the hint of the radiative precursor can be found in Ha 
emission in the zoomed profile at lower intensity level (Fig. [71 Bottom). Beyond the SNR 
boundary, continuous decreases of 11.3 fim PAH and H2 S(2) emission are seen up to the 
position around the pixel zero. The bump of the PAH emission at —4 pixel in the preshock 
area is more likely due to a background fluctuation (the bump is located outside of the 
southern boundary of Fig. [5l (b)). The profile of PAHs shows smooth decrease near the 
shock front like that of H2, but shows better agreement with that of [Ne II] beyond the 
peak to the downstream. This is consistent with the morphological characteristics of PAH 
emission seen in the line maps. Similar distribution between the 11.3 fim PAH emission 
and [Ne II] has been observed toward se veral objects. For example, in the planetary nebula. 



BD +30 3639, iMatsumoto et al.l (120081 ) detected the coexistence of PAHs and ionized gas, 
which has been attributed to a slow destruction of PAHs in the ionized region with the high 
electron density. Then, the good agreement between the 11.3 fim emission and [Ne II] in 
N49 might imply that destruction of small PAHs slowly proceeds in the shocked region (see 
more § 4.2). 

As mentioned in the previous section, the Spitzer IRAC 8 fim image is quite similar to 
the 3.3 fim PAH emission image (Fig. [5l (6)). In addition, the structure protruding from 
the bright tip at east is only seen in the 8 fim image among the IRAC images, which infers 
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that the 7 — 9 fim PAH bands possibly account for the feature. In the IRS spectra, [Ar II] 
line at 6.98 fim is the strongest emission line within the IRAC channel 4 band. However, it 
is necessary to consider other contributors such as H2 0-0 S(5), S(4) at 6.9, 8.0 /im as well as 
possible 7—9 /im PAH emissions. In order to examine their relative contributions to the IRAC 
8 /im band, we measure the intensities of H2 S(5) and [Ar II] lines. Since they are resolvable 
enough to estimate their intensities by using two Gaussian fits, we can obtain 1.76 ± 0.34 
and 3.78 ± 0.36 in unit of 10~^^ W m~^ for H2 S(5) and [Ar II], respectively. Then, the 
intensity of H2 emission together with the 7.7 fim PAH feature intensity (3.76 ± 0.52 x 10^^^ 
W m~^ in total) seems to be comparable to that of [Ar II], so that the emission seen in the 
IRAC 8 fim image is dominated by not only the [Ar II] but also the H2 lines with the 7 — 9 
/im PAH features. Therefore, the similarity between the 3.3 /im PAH emission image and 
the 8 fim image could arise from 7 — 9 fim PAH emission with H2 line emission, although 
we cannot rule out that the PAHs emitting the feature can coexist with the ionized gas. 
There is non-zero continuum level in the IRS spectrum in Figure 6 that we attributed to the 
residual of the background subtraction. The dust emission and synchrotron (or free-free) 
emission from the SNR should be negligible at this wavelength, so that even if the non- 
zero continuum is real, it should come from unrelated object. Therefore, we conclude that, 
although we cannot completely exclude the contamination from the general ISM, the good 
spatial correspondence is mainly due to the 7 — 9 fim PAH emission in the IRAC 8 fim band. 



4. Discussion 

4.1. Physical properties of PAHs in N49 

The PAH emission features are expected to reflect the local physical conditions. In 
particular, the relative strength variations among the PA H features at differ ent wavelengths 



are mainly attributed to their different charge states (e.g.. lBakes et al.ll200ll ). PAH emission 
from C-H modes such as 3.3 and 11.3 fim features becomes stronger in a neutral state, 
whereas others from C-C modes such as 6.2 and 7.7 fim features become stronger in an ionic 
state. Generally, comparisons of the intensities of 3.3 fim and 11.3 fim, both norma lized to 



that of 6.2 fim, are regarded as an indicator for the degree of ionization of PAHs (e.g.. lTielens 



20081 . and references therein). Similarly, the intensity ratios of 6.2 fim to 11.3 fim and 7.7 
fim to 11.3 fim are often used as well. In our observations, the intensities of 3.3 fim bands 
cannot be directly compared with those of other PAH bands in the IRS spectra, but from 
the IRS spectra, we calculate band ratios among 6.2, 7.7, and 11.3 fim PAH emission; 
/6.2//11.3 = 0.63 ± 0.31, /7.7//11.3 = 1-25 ± 0.26. 



Our ratios show unusually weak 6 — 8 fim PAH features relative to a 11.3 fim feature in 
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comparison with the median of the ratios from the star- forming galaxies ih.i/ Iu.?, — 3.6 in 



Smith et aLll2007D or those for Galactic diffuse emission (/7.7//11.3 — 2.0 — 3.3 in lSakon et al. 



20041 ). Such low PAH 7.7/11.3 ratios have been reported in some of elliptical galaxies, which 



is reasonably interpreted as a result of a larger fra ction of neutral PAH s due to soft radiation 
field from evolved stars (e.g., /7.7//11.3 ~ 1 — 2 in lKaneda et al.l 120081 ). Similarly, our ratios 
can indicate that the PAHs in N49 are dominantly neutral even though the environmental 
cause is different. Figure [8] is a diagram of PAH band ratios (/6.2/-^ii.3 versus /7.7//11.3) 
from various objects in literatures, which shows that these PAH band ratios are linearly 
correlated regardless of object type. Note that the PAHs are mainly neutral in the lower 
left and mainly ionized in the upper right. Interestingly, it is found that our ratios also 
follow the universal linear correlation between the two ratios. As far as we know, N49 is 
the only SNR where both IQ.2/I11.3 and /7.7//11.3 ratios are obtained. In N132D, which is 
another SNR with PAH e mission, only an upper limit of Iq.2/Iii.z ratio could have been 
derived ( iTappe et al.ll2006l ). and is marked in Figured This upper limit suggests that the 
ionization fraction of PAHs in N132D is not high, which is consistent with the theoretical 
expectation that PAHs can becorne neut ral within a short time scale in the postshock layer 
(see Figure 5 in iMicelotta et al.l l2010al ). It is known that the charge states of PAHs are 
determined by the balance between the photoionization and the electron/ion collisions. Th e 
charge state is proportional to the quantity ~ Gq^JT jrif, (for a review, see iTielensI I2OO8I ). 
where Gq is the intensity of radiation field in units of the Habing field (1.6 x 10~^ W m~^), 
T is the gas temperature, and is the electron density. The ones with small /6.2/A1.3 ratio 
(or with neutral PAHs) represent a region where Gn\/T / np is small. For a few well-studied 
photodissociation regions (PDRs), iGalliano et al.l (120081 ) have empirically interpreted the 
PAH band ratio with the above physical properties, 1^.2/ hi.s — [Go/3040 ne{cm~^)](Tgas/10^ 
K)^/2^0.53 valid in the range 400 < Go/(ne/l cm~^) (Tgas/W^ K)i/2 < 400O. 

As we will discuss in the next section, the detected PAH emission is probably either 
from the PAHs in the shocked molecular gas and/or from the PAHs in the preshock gas 
heated by radiative precursor. In dense molecular clouds, PAHs usually exist in a neutral 
(and anionic) state. Even after they experience continuous (C-type) shock waves, neutral 
PAHs predomi nate in all charged states because p hoto-ionization is not significant in these 
circumstances (IFlower fc Pineau des ForetsI l2003l ). The PAH emission could be also from 
PAHs in preshoc k gas heated by radi ative precursor. In N49, by modeling the UV and 
optical spectrum, IVancura et al.l ( 1l992l ) find that preshock densities of 20 — 940 cm~^ in the 
velocity range, 40 to 270 km s~^, are required for the denser regions with radiative emission. 
The strength of the UV precursor in radiative shocks (30 km s"^ ^ ^ 200 km s~^) is 
calculated that Gq ~ I.GtioVst, where ^,, 7 is the shock velo city per 100 km s~^ and uq is 
the preshock hydrogen nucleus density ( iMcKee et al.l 119871 and references therein). The 
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coefficient varies with the shock velocity, and it is 1.6 when Vs = 100 km s~^. Then adopting 
a preshoc k density as 150 cm~ ^ at Vg = 100 km s~^ under the requirement that Pqv"^ remains 
constant ( Vancura et aLlEood ). Go becomes 240. This UV radiation heats a column of 10^^ 
cm~^ hydrogen nu clei to ~ 5000 K a nd fully ionizes it (e.g., the case of uq = 100 cm~^ and 
Allen et allboOsK In this radiative precursor region, the ratio Goy/T/rie 



100 km s ^ in 



will be sm all 110), so that the PAHs will be neutral. If we extrapolate the empirical 
relation of iGalliano et al.l (120081 ) . we obtain /6.2/A1.3 — 0.53 which is consistent with the 
observed ratio. 



Meanwhile, the observed band ratio can be compared with theoretical studies by lDraine fc Li 



(I2OOII ). which show how the charge state and size of PAHs affect relative strengths of three 
band emission, 6.2 /im/7.7 fim versus 11.3 fim/7.7 fim (see Figure 16 in their paper). Our 
observed fluxes give lQ.2/h.7 — 0.50 ± 0.27 and /n. 3/17.7 — 0.80 ± 0.15 which are located 
near the line for neutral PAHs. More interestingly, our high I&.2/h.7 ratio indicates that 
small PAHs are dominant in N49 although the error is large. The existence of small PAHs 
can also be verified by /3.3//11.3 ratio since larger PAHs produce rel atively a strong 11. 3 
ixm feature while the 3.3 /im emission is enhanced by small PAHs (e.g., ISchutte et al.lll993l ). 
It is not simple to directly derive the /3.3//11.3 ratio from our data because the 3.3 /xm and 
11.3 /im fluxes are observed by different telescopes and positions. However, we make a rough 
estimation to obtain both the surface brightness at the peak and the integrated intensity 
with the full aperture (50.4") along the IRS slit from the AKARI 3.3 /im line map (Fig. 
E]) and the Spitzer spectrum (see the slit position in Fig. E]). The /3.3//11.3 ratios of the 
peak brightness and the total intensity are measured as ~ 0.43 and ~ 0.30, respectively. 
These ratios are in a high end of the observed range (in general, 0.2 — 0.3), which can 
be accounted for by small PAHs unless the radiation is hard enough to excite large PAHs 
( Mori et al.llin prep.l ). For 100 km s~^ radiative shocks, the UV radiation from the shock 



is mostly in Lya photons (IMcKee et al.l 119871 ). so it might be difficult to directly compare 
the shock radiation to the radiation from a central source. However, no detectable PAH 
feature at longer wavelengths easily produced by large PAHs (e.g., 16.2 or 17.4 /im features) 
is more likely to support a sm all contribution of large PAHs to the 3.3 /xm feature (see Fig- 
ure 9 of IWilliams et al.ll2006l ). Therefore, the observed 1^.2/17.7 and /3.3//11.3 ratios indicate 
the presence of small PAHs, which is consistent with the PAH formation by fragmentation 
from larger carbonaceous grains but not with the preferential destruction of small PAHs in 
shocked gas. 
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4.2. Origin of PAH emission 



PAH emission has been detected toward few shock-associated regions. Due to the lim- 
ited observational evidence, it is poorly understood how shocks play a role in PAH pro- 
cessing. Generally, there are two sides to SNR shocks in terms of PA H processing, i.e ., 



shocks can produce PAH molecules by shattering larger dust grains (e.g.. lJones et al.lll996l ). 
but shocks can also de st roy them by col l isions w ith energetic particles in shocked gas (e.g., 
Micelotta et al.ll2010al ). iMicelotta et al.l (j2010bl ) show that the destruction mechanism in a 



hot gas (T > 3 X 10^ K) is dominated by electron collision for small/medium size PAHs 
(50 — 200 C-atom). They describe a PAH lifetime in a hot gas as tq = Nq/ Jrin/c^ where 
Ac is the initial PAH size, J is the rate constant for electrons, nuclear and electronic in- 
teractions, and TT-H/e is the hydrogen/electron density. At a fixed PAH size, a PAH lifetime 
can be determined by the gas density and temperature. In a rarefied gas with fast shocks 
{vs ^ 300 km s~^), no PAH formation is expected, and most PAHs are completely destroyed 
soon after being swept up by shocks. In particular, smaller PAHs are eroded more easily for 
T < 10^ K. Nevertheless, if we see the emission from PAHs just swept up, the PAH emis- 
sion should be confined to the shock front. In a dense medium, however, as shocks become 
slower {vs ^ 1 00 — 150 km s~^), PA Hs might survive and/or be formed in those shocked 
environments (IMicelotta et al.ll2010al l. Then, the detection of PAH emission in the shocked 
regions requires either efficient formation of PAHs or shielded condition against complete 
destruction. Currently, observational difference between newly formed and survived PAHs 
is not obvious, but it is supposed that the newly formed PAHs undergo the same processing 
as the survived ones soon after. Meanwhile, PAH emission could be also produced outside 
the shocked region by radiative precursors of shocks. Hence, it is necessary to separately 
examine the dominant mechanism for PAH emission depending on the local shock condition. 

The complex structures of N49 imply that physical conditions inside the remnant should 
significantly vary in localities. This remnant can be largely divided up into three kinds of 
environments related to the characteristics of the medium from which different emission 
arises. The hot gas with X-ray emission is generated by the fast shock propagating in the 
diffuse intercloud medium (preshock density uq = 0.9 cm~^, Ug = 27—2300 cm~^, T > 7x 10^ 
K, and Vg ~ 730 km s~^ in lPark et al.ll2003l ). The dense ambient medium produces the bright 



optical filame nts (i.e., ionic line em ission) by shocks of v < 140 km s~^ (tt-q = 20 — 940 cm~^, 
T ~ 10^ K in lVancura et al.lll992l ). The electron densities of the filaments var y from Up ^ 70 
to 515 cm"'^, and the dens est regions have the densities of 1000 — 1800 cm~^ (IBilikova et al. 



20071 : IVancura et al.lll992l ). Lastly, the molecular line emission such as H2 lines can arise 
from the dense clumps where slo w (nondissociative) sho cks are impinging on (in general. 



f^o ^ 30 cm ^, Vs < 50 km s ^ in lDraine fc McKeelll993l ). In N49, both shocked molecular 



and ionic gases interestingly show analogous distributions, and the PAH emission is likely 
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to be associated with both gases. This is not surprising; in the s ense that N49 is known to 
be interacting with a nearby molecular cloud (IBanas et al.lll997l ). The coexistence of both 
molecular and ionic shocks in a single SNR is occasionally seen in SNRs interacting with 
molecular clouds, since the shocked dense clumps emitting H2 transitional line emission are 
immersed in a less dense rn edium that emits various ionic line emission after swept up by 
shocks (e.g.. IChevalierlll999l ). 



In the hot diffuse medium of N49 (rig = 27 cm ^, T = 7 x 10^ K), the lifetime of 
PAHs with Nc = 200 is ~ 1 month by adopting the analyt ical fits to the rate constant 
(J ~ 2.7 X 10"^ cm^ s"^ Table 2 in iMicelotta et al.l boiObl ). Even for large-sized PAHs 
or PAH clusters, the li fetime is not long enou gh to be detected (e.g., ~ 20 months for 
Nc = 1000, Fig. 8 in IMicelotta et al.l l2010b( ). To explain the observed PAH emission, 
protective environments in N49 are necessary, which have already b een noticed in optical 



spectroscopic observations (IVancura et al.lll992 



Bilikova et al. 



2007). By using the echelle 



observations with the HST image (see Fig. 3 in lBilikova et al.l 120071 ). some broad emission 
features in the echellegrams reveal a "head-tail" structure, of which the bright emission at the 
head has a smaller radial velocity offset from the systemic velocity than the fainter emission 
in the tail. This type of the structure could originate from shocks encountering the dense 
gas. While a shock is propagating through the preshock gas in the inner part of the dense 
bright regions, the shock velocity drops so that the gas deep inside the dense filament can 
avoid to be swept up by shocks, staying unshocked. The observed Ha recombination lines 
at the systemic velocity in the optical spectrum support this. 

Based on the above discussion, we propose a schematic model of the PAH emitting 
condition in N49 in Figure [H As a shock is retarded by a dense gas (Region II and shocked 
Region HI), the temperature of the postshock gas can substantially decrease (Tgas oc v^). The 
PAH lifetime under these circ umstances can significa n tly increase because t he rate constant 
rapidly declines below 10^ K (IMicelotta et al.ll2010bl ). IVancura et al.l (119921 ) have measured 
electron densities and electron temperatures for the bright optical filament of N49 as ~ 1 x 10'^ 
cm~^ and ~ 1 x 10^ K, respectively. These quantities lead to the lifetime of 1.2 x 10^ yr 
for 200 C-atom PAHs (c.f., J ~ 5.2 x 10~^^ cm^ s~^), which is much longer than the age of 
the SNR. PAHs in the central region (unshocked Region HI) of the molecular clumps (and 
the dense filaments) can be heated by UV photons that mainly originate from the radiative 
shell, more precisely in this case, the knotty optical filaments around the clumps. In the case 
of the PAH emission observed outside the shock front (seen eastward at the tip in Fig. [3]), 
it is most likely that the preexisting PAH molecules are heated by the radiative precursors 
and produce the emission features. When the shock velocity is greater than 80 km s~ ^, the 



ionizing UV fluxes produced by the precursor become effective (IShuU fc McKee 



bulk range of the shock velocity in N49 is measured as < 140 km s ^ (IVancura et al.lll992l ) 



1979). The 
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and we indeed see the hint of the radiative precursor in the zoomed Ha hne profile (Fig. [7]). 
In summary, PAH emission can be associated with shocked ionic gas (Region II) or (shocked) 
molecular gas (Region HI) when the shock is sufficiently retarded with existence of heating 
source such as UV photons (Fig. Also, as a shock is propagating and retarded inside the 
dense region, shocked ionic and molecular gas can hierarchically exist together with PAH 
emission. Then, the morphological correlation of the PAH emission to other emission would 
depend on the shock velocity and the preshock density. 

The absence of PAH emission at the peak of the H2 emission is interesting. In Figure 
[5l the region with the brightest H2 emission seems to coincide with one of the bright Ha 
filaments, and it is found that the filament particularly shows the most extreme difference 
among Ha, [S II], and [O HI] emissions in optical; the [O HI] emission is absent while the [S 



II] em ission is relatively stronger than Ha (see Fig. 1 i n iBilikova et al.l 120071 : IVancura et al. 



1992 ). According to shock models in Vancura et al. ( 19921 ). [O HI] emission is generally 



produced by faster shocks and cannot be produced for shocks with velocities < 80 km 
s~^. Meanwhile, as the shock velocity decreases, the intensities of [S II] show an increasing 
tendency with respect to that of H/3 in the model calculation. In this context, it is likely that 
relatively slow shocks (< 80 km s~^) have reached the filament and have resulted in these 
differences in the optical spectra. If so, the faint PAH emission can be ex plained, because 
the s hocks with Vs < 80 km s~^ cannot produce sufficient UV radiation f lShull fc McKee 



I979I ) at the outer layer of the filament which is necessary to heat PAH molecules. For 



H2 emission, however, hydrogen molecules can be collisionally excited by shocks with those 
velocity ranges, and produce bright emission. For Vg > 25 — 30 km s~^, the li ne intensity is 
sensitive to the preshock density, not to the shock velocity (IBurton et al.lll992l ). Hence, even 
if PAHs can survive, they might not generate detectable emission due to the lack of heating 
sources, while it is possible to produce bright H2 emission in the same environment. This 
may also explain the IRS spectra of the Galactic SNRs interacting w ith molecular clou ds, 
where H2 emission is prominent but the PAH emission is absent (e.g.. iHewitt et al.ll2009l ). 



The detection of 3.3 /xm PAH band feature from N49 is surprising in a sei ise that only 
small PAHs with < 100 C-atoms or a size < 6 A can notably produce the feature ( Draine fc Li 



20071 ) while small PAHs can be easily destroyed by strong shocks. Even elliptical galaxies 



showing low PAH 7.7/11.3 fim ratios (i.e., mainly neutral PAHs) seem to have no significant 
3.3 fim PAH emission due to the lack of small PAHs (IKaneda et al.l 120071 ). Meanwhile, 
the com parison between our PAH band ratios {I&.2/h.7 and /3.3//11.3) and the numerical 



studies ( jPraine fc Lill2001t iMori et al. 



in prep. I ) support that small PAHs are dominant in 
N49. In addition, there is no indication of the 15 — 20 /im feature in the IRS spectrum by 



Williams et al.l ( l2006l ) (see their Fig. 9), which is a signature of large PAHs. In the case 
of N132D, however, the relatively strong 15 — 20 /im hump compared to 6 — 11.3 /im PAH 
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features suggests that the large PAHs are domin ant, which was int erpreted as an evidence 
for the survival of large PAHs behind fast shock flTappe et al.ll2006l ). Thus, PAH processing 
in N49 and N132D is likely to differ from each other. This difference in PAH sizes can be 
explained by the different envi ronments and ev olutionary stages of the two SNRs: N49 is a 
middle-aged SNR (~ 6600 yr, iPark et al.l 120031 ) and is interacting with a nearby molecular 
cloud with a radi ativ e shock {vg — 100 km s~^), while N132D is relatively young (~ 2500 
996) and still has a fast shock in less dense environment (v^ ^ 800 km 



yr 



Morse et al. 



Morse et al. 



19961 ) ■ In a dense gas with a moderately low temperature 3 x 10^ K), 



PAHs would survive regardless of the size. This is because the nuclear interaction with 
helium is the dominant process of destructi on in that situation, wh ich does not significantly 
depend on the size of PAHs (see Fig. 5 in lMicelotta et al.ll2010bl ). Hence, small PAHs are 
NOT preferentially destroyed in N49, which can reasonably explain our detection of the 3.3 
/im PAH feature in the postshock region with the relatively low temperature. This might 
imply that a wide range of PAH properties can be found in SNRs. 

Based on our observations, it might be difficult to understand in detail, for example, 
what is a dominant heating mechanism for emitting PAH features in the SNR, and whether 
the PAHs are newly formed or not. Nevertheless, it seems obvious that one of the essential 
conditions for the survival/(re)formation of PAHs is the existence of dense ambient medium 
around a SNR which directly influences the evolution of shocks such as the shock velocity. In 
addition to this, the sufficient heating source such as UV radiation is necessary for survived 
PAHs to produce observable PAH emission from a SNR. 



5. Summary 

We have carried out an infrared (IR) spectroscopic study on the supernova remnant, 
N49 in the Large Magellanic Cloud by using AKARI IRC NG observations which cover most 
of the bright eastern regions including distinct filaments. Since the observations have been 
performed as a coarse spectral mapping, we are able to make spectral line maps and have 
compared the distribution of different emission features. In the AKARI IRC spectra (2.5 — 
5 /um), we detect 3.3 /im polycyclic aromatic hydrocarbon (PAH) band features with several 
strong hydrogen recombination lines and moderate H2 molecular lines. To our knowledge, 
this is the first time that we observe the presence of the 3.3 fim PAH feature related to a 
SNR. Our main results are summarized below: 

1. The 3.3 fim PAH feature in the spectra are clearly distinguished from other shocked 
ionic/molecular lines in terms of the line width and the intensity variation with position. In 
the AKARI line maps, the distribution of the 3.3 /im PAH emission shows overall association 
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with those of other emissions such as Bra and H2 1-0 0(3), which indicates that the PAH 
emission indeed originates from the SNR and is associated with both ionic and molecular 
emission. In addition, there are morphological dissimilarities among the line maps in a 
local scale reflecting the different physical conditions such as shock velocity, and preshock 
density in each region. The morphological characteristics are also clarified and confirmed by 
comparison to archival Spitzer 8 /im, HST Ha, and H2 2.12 /im images obtained at AAT. 

2. Overall distribution of the 3.3 /im PAH emission is more similar to that of the 
H2 emission in a sense that the H2 emission is spatially confined but the Bra emission is 
extended over a large area. In addition, the 3.3 /im PAH emission is extended beyond the 
eastern shock boundary with the intensity comparable to those in the inner regions. These 
indicate that the PAH emission is possibly associated with the interacting molecular cloud, 
the center of which nearly coincides with the bright PAH emission region. Meanwhile, the 
peak position of the PAH emission is not matched with that of the H2 emission but rather 
close to that of the Bra emission. In addition, the PAH emission is relatively faint at the 
peak of the H2 emission. 

3. We find signatures of other PAH features, C-C stretching modes at 6.2 and 7.7 
/im and C-H out-of-plain bending mode at 11.3 yum in the archival Spitzer IRS SL spectra 
(5.2 — 14.5 /im). We derive the band ratios of PAHs, IQ.2/I11.3 and /7.7//11.3 from the IRS 
spectra (/6.2/A1.3 = 0.63 ± 0.31, /7.7//11.3 = 1.25 ± 0.26), which implies that the PAHs in 
N49 are dominantly neutral. This is consistent with the theoretical expectation for shocked 
PAHs in dense molecular clouds. It is found that the ratios follow the universal linear 
correlation between the two ratios in literatures. We try to asso ciate the /fi.2/-^i i ..s rat io with 



the physical quantity Goy/T/rie by using the empirical relation in lGalliano et al.l ( 20081 ) under 
the assumption that the shock radiation is the dominant heating mechanism of PAHs in a 
very dense clump. In addition, the relatively high lQ.2/h.7 and /3.3//11.3 ratios (0.50 ± 0.27 
and ~ .36, respectively) i ndicate the existenc e of small PAHs according to the numerical 



studies (jPraine fc Lill200ll : iMori et al.llin prep. I ). These results are consistent with the PAH 



formation by fragmentation from larger carbonaceous grains but not with the preferential 
destruction of small PAHs in shocked gas. 

4. The morphological features of the PAH emission can be attributed to the different 
mechanisms of SNR shocks in terms of PAH processing. For the PAH emission associated 
with either shocked H2 gas or (shocked) ionic gas, or both, PAHs must exist in dense gas 
where the shocks have been sufficiently retarded (even terminated) to avoid the complete 
destruction. Depending on the shock velocity and the preshock density, the PAH emission 
can be associated with either the ionic gas or the molecular gas. Although PAHs can survive 
a slow shock, detectable PAH emission may not arise due to the lack of UV radiation in a 
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certain condition. For the PAH emission outside the SNR, the radiative precursor could be 
responsible for the excitation. For PAHs to exist and radiate in SNRs, an ambient dense 
medium and a sufficient heating source around the medium are most likely to be required. 

This work is based on observations with AKARI, a JAXA project with the participation 
of ESA. We wish to thank all the members of the AKARI project. We would like to thank 
John Dickel for providing the AAT data and E. R. Micelotta for helpful discussions about 
PAH processing in a shocked region. We are also grateful to J.-H. Shinn for useful interactions 
about the calculations of H2 level population. This research was supported by Basic Science 
Research Program through the National Research Foundation of Korea (NRF) funded by 
the Ministry of Education, Science and Technology (NRF-20 11-0007223). 
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Fig. 1— Positions of the IRC NG/Ns slits superposed on the AKARI N3 image of N49. 
Thirteen on-positions and two background positions are observed (black). We also mark 
where final averaged spectra are extracted (white). The five final spectra are referred to as 
"PI" to "P5" from east to west, respectively (see Fig. [2]). For comparison, the contours 
showing 6-cm radio continuum ( green contours) a nd CO emission (red contours) are super- 



posed. The radio data are from iDickel fc Milne I (119981 ): the levels are at 2, 4, 6 , 8, and 10 
mJy beam~^. The CO contours are taken from Figure 5 of iBanas et al.l (119971 ): the levels 
are at 1.0, 1.3, 1.7, 2.0, 2.3, and 2.7 K km s"^. 
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Fig. 2. — Final spectra together with the background spectrum. For convenience, the spectra 
are shifted by 1, 2, 3, and 4 xlO~^ W /im~^ sr~^ for P2 to P5, respectively. Background 
spectrum is shown in the bottom (shifted by -1 xlO~^ W m~^ nm~^ sr~^). Noticeable lines 
are labeled at their vacuum wavelengths. 
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Fig. 3.— Ha 1-0 0(3), PAH, and Bra line maps of N49 by using AKARI IRC spec- 
tra. For comparison , contours (black) are taken from an Ha image of the HST WFPC2 
(IBilikova et al.l 120071 ): these contours are 10, 30, 50, and 70 % of the peak. The black cross 
marks the peak of PAH emission at (aaooo, ^2000) = (05'*26™05.14^ -66°05'03"). The white 
crosses denote two clumps seen in the PAH line map. The scale bar represents 10" or 2.4 pc 
at 50 kpc. The color bars are given in units of W m~^ arcsec"^. North is up, and east is to 
the left. 
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Fig. 4. — Variation of line intensities and H2 column densities with respect to the different 
positions (from PI to P4). The estimation of the column densities are described in § 3.2. P5 
spectrum is excluded due to non-detection of any H2 lines. 




Fig. 5.— Images of N49 in f a) 2.12 jPickel et al.l ll995h. (b) IRAC 8.0 /im, (c) and (d) 



smoothed and orig inal Ha (iBilikova et al.ll2nn7[ K (c) Ha image is convolved to the spatial 
resolution of AKARI N3 image (FWHM ~ 4.0"). Contours overlaid in all images are from 
the PAH 3.3 /im map of Figure |3l of which levels are at 6, 7, 8, and 9 xlO~^^ W m~^ 
arcsec"^. North is up, and east is to the left. Positions of the Spitzer IRS SL resolution 
slits are marked in (b) and (d). 
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Fig. 6. — (a): IRS SL spectra of N49 with Spitzer IRAC spectral response curves. "Nod 1" 
position spectra are used only. The [Ne II] line at 12.8 /im is truncated for clarity. Dominant 
emission lines are labeled. (6) — (d): Profiles of the PAH band emission features at 6.2, 11.3, 
and 7 — 9 extracted from the IRS spectra. The extracted ranges are overplotted in (a). 
Detected PAH bands with H2 emission lines are also marked. Error bars at the left bottom 
of each panel represent the typical error for the spectra. 
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Fig. 7.— Top: IRS SL line profiles of PAH A 11.3 /im (black solid), Hs 0-0 S(2) A 12.3 
/iin (dotted), and [Ne II] A 12.8 fim (dash dot). For comparison, an one-dimensional cut of 
Ha emission along the IRS slits (gray solid) is extracted from the HST image smoothed to 
the pixel scale of the IRS. Bottom: The same profile of Ha, but zoomed in for the variation 
at lower levels. The dashed line represents the SNR boundary (the probable location of the 
shock front). One pixel corresponds to 1.8", and the coordinate of the position for the zero 
pixel is (a, 6) = (05''26"*05", -66°05'28"; J2000.0). The negative and the positive correspond 
to south and north, respectively. 
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Fig. 8. — PA.H band ratios of SNR N49 as well as various objects from liter at ures 



fIVermeii et al.ll2002l : ISloan et al.ll2007l : iGalliano et al.ll2008l : iBernard-Salas et al.ll2009[ ). The 
ratio in N49 is derived from the Spitzer IRS spectr um (j^ 3.3). The up per limit of h.2/Iii.z 
from the SNR N132D is marked as a dashed line (ITappe et al.l 120061 ) . Different symbols 
denote different types of objects. Few objects without any symbols are either Herbig Ae/Be 
or T Tauri stars or a Magellanic object. 
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Fig. 9. — Schematic diagram of local medium structures in N49. There are largely three 
kinds of medium; the hot interior, the radiative shell, and the molecular clumps emitting 
X-ray, ionic line emission, and molecular emission, respetively. Darker color indicates the 
higher survivability of PAH molecules. In the case of a dense molecular clump, its inner 
layer can remain unshocked. UV photons mainly originate from the radiative shell, more 
precisely, neighboring optical filaments around the molecular clump in N49. Survived PAHs 
associated either the radiative shell or the molecular clump can be heated by UV photons 
and produce detectable emission features. 
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Table 1. Summary of Observations 



Obs. ID 


Sub. ID 


Program^ 


Date 




Field center 
R.A. (J2000) Dec. (J2000) 


Note^ 


1420813 


1 


ISMGN 


2008 


Jun 


12 


05:26:12.1 


-66 


05:26.8 


BG 


1420813 


2 


ISMGN 


2008 


Jun 


23 


05:26:12.2 


-66 


05:27.3 


BG 


1420925 


1 


ISMGN 


2008 


Jul 


10 


05:26:04.9 


-66 


05:12.2 


P2 


1420926 


1 


ISMGN 


ZUUo 


Till 
JUl 


08 


UO.ZO.UO.O 


-00 


uo.uu.y 


P3 


1420926 


2 


ISMGN 


2008 


Jul 


08 


05:26:04.0 


-66 


04:56.7 


P3 


1421925 


1 


ISMGN 


2008 


Jul 


10 


05:26:05.7 


-66 


05:16.6 


PI 


1421925 


2 


T TV IT TV T 

ISMGN 


2008 


Jul 


10 


05:26:05.5 


-66 


05:14.7 


PI 


1421926 


1 


ISMGN 


2008 


Dec 


20 


05:26:02.6 


-66 


05:04.5 


P4 


1900451 


1 


P3LMC 


2008 


Jun 


25 


05:26:04.6 


-66 


05:13.8 


P2 


1900452 


1 


P3LMC 


2008 


Dec 


25 


05:26:03.4 


-66 


05:12.1 


P3 


1910233 


1 


LMCNG 


2009 


Dec 


11 


05:25:58.9 


-66 


05:16.1 


N 


1910233 


2 


LMCNG 


2009 


Dec 


12 


05:25:56.8 


-66 


06:15.4 


N 


1910233 


3 


LMCNG 


2009 


Dec 


12 


05:26:00.6 


-66 


05:08.7 


P5 


1910234 


1 


LMCNG 


2009 


Dec 


13 


05:26:00.3 


-66 


05:13.4 


P5 


1910234 


2 


LMCNG 


2009 


Dec 


13 


05:26:00.5 


-66 


05:18.9 


P5 


1910234 


3 


LMCNG 


2009 


Dec 


15 


05:26:00.6 


-66 


05:11.1 


P5 



Note. — 

"^Name of the AKARI mission programs: "ISM in our Galaxy and nearby galaxies" 
(ISMGN, PI: H. Kaneda), "Observations of the Magellanic Cloud" (P3LMC, PI: T. On- 
aka), and "Near infrared spectroscopic observations of red objects in the Large Magellanic 
Cloud" (LMCNG, PI: T. Onaka). 

^Classification of the final spectra after averaging separate spectra extracted from over- 
lapping region. Spectra with the same note are used for one final spectrum. BG: back- 
ground spectrum, P1-P5: source spectra extracted from five different regions, respectively 
(see Fig. [1]). N: spectrum not used for the final spectra owing to either a single observation 
under a poor detector condition or a wrong pointing. 
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Table 2. Detected lines and their intensities from AKARI IRC Data for N49 



Intensity (10 ^ W m ^ sr 



T irip IT) 


^ref 


PI 


P9 


P'^ 


P4 


P^ 


Br/3 


2.626 


<0.34 


1.58±0.30 


2.73±0.27 


2.47±0.53 


0.77±0.34 


H2 1-0 0(3) 


2.803 


0.27±0.15^ 


0.31±0.25 


1.38±0.23 


0.85±0.40 


<0.47 


H2 1-0 0(5) 


3.235 


0.32±0.15^ 


0.75±0.38 


0.58±0.12^ 


0.42±0.24*^ 


<0.57 


Pf5+PAH^ 


3.297 


1.72±0.31'^ 


2.16±0.51 


1.42±0.25'i 


0.79±0.37'i 


<0.73 


PAH<^ 


3.300 


1.67±0.35 


1.87±0.53 


1.04±0.28 


0.39±0.45 


— 


PAH^ 


3.400 


— 


0.61±0.36f 


0.80±0.34f 


— 


— 


Pf7 


3.741 


<0.24 


0.22±0.09^ 


0.51±0.13 


<0.44 


<0.45 


Bra 


4.052 


0.50±0.16 


2.94±0.15 


3.83±0.13 


4.03±0.25 


1.52±0.17 


H2 0-0 S(ll) 


4.181 


<0.09 


0.23±0.08^ 


0.50±0.11 


0.42±0.14^ 


<0.48 


H2 0-0 S(10) 


4.410 


<0.20 


0.20±0.09'^ 


0.35±0.14 


<0.51 


<0.29 


Pf/3 


4.654 


<0.57 


0.84±0.11^ 


0.79±0.11^ 


0.88±0.17^ 


0.64±0.24 


H2 0-0 S(9) 


4.695 


<0.44 


0.26±0.20'^ 


1.29±0.23'^ 


1.17±0.35'^ 


<0.55 


[Fe II] 


4.889 


<0.35 


0.40±0.16^ 


0.53±0.23 


<0.48 


0.49±0.18*^ 










Ratio 








Case B^ 


PI 


P2 


P3 


P4 


P5 


Br/3/Bra 


0.60 


<0.68 


0.54±0.11 


0.71±0.07 


0.61±0.14 


0.51±0.23 


Pf^^^ /Bra 


0.10 


3.44±1.26 


0.73±0.18 


0.37±0.07 


0.20±0.09 


<0.48 


Pf7/Bra 


0.135 


<0.48 


0.07±0.03 


0.13±0.03 


<0.11 


<0.30 


Pf/3/Bra 


0.20 


<0.57 


0.29±0.04 


0.21±0.03 


0.22±0.04 


0.42±0.16 



Note. — Intensities with errors of detected lines are given. Upper limits are at the 
2a level. 

'^Vacuum wavelength of a line in /xm. 

^FWHM fixed to the instrumental width (~ 0.03 yum) for a gaussian fit. 

■^Total intensity of the feature at 3.3 /im, which can be mainly contributed by either 
3.3 yum PAH or Pf6 emission. 
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'^Central wavelength fixed for a gaussian fit. 

*^Estimation of PAH intensity by subtracting the total intensity of the 3.3 /^m feature 
by the scaled Bra intensity (i.e. 10% of the Bra intensity). 

^Intensity of a minor PAH feature at 3.4 /xm. Because of the curvature around the 
feature, the baseline is fitted by using a second order polynomial. Due to the ambiguous 
identification, the measured intensities can be regarded as an upper limit of the 3.4 
/im feature intensity. 



^Theoretical recombination line ratio in the case of "On-the-spot" approximation (so 
called "Case B" ) with Ue ~ 10^ - 10^ cm'^ and Te ~ 5 - 30 x 10^ K jHummer fc Storev 
1987h . 



